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Abstract. We calculate multicomponent line-driven wind models of stars at extremely low metallicity suitable for massive first 
generation stars. For most of the models we find that the multicomponent wind nature is not important for either wind dynamics 
or for wind temperature stratification. However, for stars with the lowest metallicities we find that multicomponent effects 
influence the wind structure. These effects range from pure heating to possible fallback of the nonabsorbing wind component. 
We present a simple formula for the calculation of metallicity for which the multicomponent effects become important. We 
show that the importance of the multicomponent nature of winds of low metallicity stars is characterised not only by the low 
density of driving ions, but also by lower mass-loss rate. 
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1. Introduction 

First generation stars are a textbook example of the importance 
of initial metallicity for stellar structure and evolution. During 
the gravitational collapse of extremely-low metallicity clouds 
the formation of massive stars was possibly much more favored 
than in the present time (cf. Bromm et al. 2002, Nakamura & 
Umemura 2002 1 and the evolution of such massive stars was 
strikingly different from standard stellar evolution (Siess et al. 
2002 1. Similarly, due to their enormous luminosity and minute 
abundance of elements heavier than helium, their line-driven 
winds substantially differed from those of present day stars. 

The first models of line-driven stellar winds suitable 
for very massive first generation stars were calculated by 
Kudritzki (2002). He concluded that these stars must be near 
to the Eddington limit in order to possess line-driven winds. 
However, in the case of stars at such extremely low metallicity 
other effects may become important. It is well known that ra- 
diative driven stellar winds have a multicomponent nature (cf. 



Castor et al. 1976), the reason for which is relatively straight- 
forward. The radiative force accelerating a stellar wind is dis- 
tributed unevenly over individual ions. Whereas minor absorb- 
ing elements, such as C, N, O or Fe, obtain momentum from 
the stellar radiation field, hydrogen and and helium are only 
marginally accelerated by the radiation. However, these nonab- 
sorbing components are accelerated by friction with the absorb- 



ing components. This momentum transfer between low-density 
absorbing and high-density nonabsorbing components is spe- 
cially important for low-density stellar winds (Springmann & 
Pauldrach [T992l Krticka & Kubat 2001 hereafter KKII). These 
stellar winds can be heated by the friction between components 
and, for extremely low wind density, even hydrogen fallback 
may occur. Low density stellar winds may also be subject to so 
called runaway instability (see Owocki & Puis 2002 Krticka 
& Kubat 2002). This instability occurs when the velocity dif- 
ference between absorbing and nonabsorbing components is 
comparable to the averaged sound speed (Krticka & Kubat 
2001, Eq. (45)). Finally, low-density stellar winds are heated 
by the so-called Gayley-Owocki heating (Gayley & Owocki 
1994, hereafter GO). This heating/cooling is caused by the de- 
pendence of the radiative force on the velocity via the Doppler 
effect. 

Because these multicomponent effects usually occur 
mainly in the outer part of the wind downstream from the crit- 
ical point they do not affect the mass-loss rate but they influ- 
ence the wind temperature and outflow velocity. In this paper 
we discuss these multicomponent effects for stars at extremely 
low metallicity. 

2. Domain of importance of multicomponent 
effects 
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In order to assess the importance of multicomponent ef- 
fects we use a simple formula for the velocity difference be- 
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tween absorbing and nonabsorbing components. Similarly to 
Springmann & Pauldrach (IT9921 and Owocki & Puis d2002l 
we start from the equation of motion of nonabsorbing compo- 
nent for a two-component wind. Neglecting the gravitational 
acceleration, electrical force and gas pressure term, this equa- 
tion takes the form (see KKII) 



dv r 



2„2 



dr 



m p m; kT 



lnAG(xpi), 



(1) 



where v rp is the velocity of the nonabsorbing component. For 
the calculation of frictional acceleration on the right-hand side 
of Eq.([0 we assume that the temperatures of all components 
are nearly equal (T p « T; « T), and p\ is the mass density 
of the absorbing component, m p , m u q p and q t are the parti- 
cle masses and charges of nonabsorbing and absorbing compo- 
nents, and In A is the Coulomb logarithm. The Chandrasekhar 
function G(x pl ) is defined in terms of the error function erf (x) 
(cf. Burgers [Tic39l 



G(x pi ) 



erf(xpi) 



exp i 



(2) 
(3) 



The Chandrasekhar function can be approximated for x p i < 1 
by 

GW - f% (4) 

Using a /3-velocity law v rp = foo(l — 4r) (with (3=1, which 
is near to the mean observed value, cf. Puis et al. 1996 1 for the 
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non-absorbing component we can approximate 
The continuity equation can be used to calculate the density of 
absorbing ions 
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where and are the terminal velocity and stellar radius, 
2)i is the mass density ratio of absorbing and nonabsorbing 
ions in the solar photosphere (2)j = 0.0127, this value cor- 
responds to the solar ratio of sum of densities of C, N, O, Fe to 
the density of bulk plasma), Z/Zq is the metallicity (number 
density of absorbing ions relative to hydrogen) in the stellar at- 
mosphere relative to the solar value, and SOT is mass loss rate. 
Solving the momentum equation Q for the velocity difference 
we obtain 



VooR* i\pKm p m,ikT 
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(6) 



Apparently, as the wind is accelerated the velocity difference 
increases. Thus, multicomponent effects are important mainly 
in the outer part of the wind. If the wind temperature is constant 
then the maximal velocity difference is attained for maximal 
radial velocity, ie. for v rp = Voo. This enables us to calculate 
the maximal velocity difference between wind components for 
each model star. 
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Fig. 1. The domain of importance of multicomponent effects. 
Each model star is represented by a point in this plot. The 
straight line (calculated using Eq.(|6j, where we inserted v rp = 
Doo) divides regions where multicomponent effects are negligi- 
ble from where they are important. Winds of stars above this 
line have a low velocity difference between wind components 
whereas for stars below the line multicomponent effects are im- 
portant. Different symbols denote different stellar luminosity 
(given in solar units). 



The multicomponent effects are found to be important 
when the nondimensional velocity difference is typically 



== > 0.1 



(KKII). In this case frictional heating influences the wind tem- 
perature. Note that there is some kind of feedback because 
the higher temperature enhances the multicomponent effects. 
In the plot of metallicity versus the quantity v^R^/Wl (see 
Fig^ we shall use a straight line to divide regions where 
multicomponent effects are important/unimportant. Stars above 
this line have a low relative velocity difference between wind 
components (see Eq.(|6j) whereas stars below this line have 
higher relative velocity difference and thus multicomponent ef- 
fects are important for these stars. The dividing line (for which 
Xpi = 0.1) was plotted with assumed average wind temperature 
T = 25 000 K, ionic charges q p = 0.85e, q^ = 4.0e (where e 
is the electronic charge) and the mean mass of absorbing ions 
corresponding to carbon. 



We added each star for which Kudritzki ( 2002 ) calculated 
a wind model into Fig^ We used the same wind parameters 
(ie. wind temperature, ionic charges and mass of absorbing 
ions) as in the previous paragraph for the plot of dividing line. 
There are several stars for which the multicomponent effects 
are important. For stars with highest luminosity only the mod- 
els with lowest metallicity may suffer from the multicompo- 
nent effect. The lower the luminosity, the lower the mass loss 
rate and the higher is the metallicity for which multicompo- 
nent effect becomes important. Moreover, for stars with given 
basic parameters (ie. mass, radii and effective temperature) we 
are able to find a value of metallicity for which the multicom- 
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ponent effects influence the wind structure. This conclusion is 
especially important for low metallicity stars. 

For normal solar-type ion abundance with (Z/Z Q ) w 1 Eq. 
(|6j indicates that ion runaway should occur only for very low 
mass loss rates, i.e. 97t ~ 10 -11 M Q /yr, and so is relevant only 
for relatively low luminosity stars. But for stars with very low 
metallicity, ion runaway could become important for quite lu- 
minous stars. Consider, for example, the canonical CAK mass 
loss scaling (see Pauldrach et al. 1986 1 
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where T is the Eddington parameter, and Q sa 2000 Z/Zq 
(Gay l ey [T995V Because a < 1 (Abbott H982I Puis et al. 
2000 1 this generally means that lower metallicity enables larger 
multicomponent effects not only due to the lower abundance 
of the driving ions but also due to the lower mass-loss rate 
(see Eq.(|(j)). The mass-loss rate formula ([7) can be rewrit- 
ten using the scaled quantities M_ n = 97t/(lO _11 OJt0/yr), 
M 2 = an/10 2 OTt Q , L 6 = L/lO 6 i , and T w 0.26L 6 /M 2 . 
For the canonical case of a = 2/3 and T = 0.5, this yields 
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Similarly, Eq.0 implies that the maximal velocity difference 
scales as 
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where v s = u oo /(10 8 cms _1 ), R 12 = i?*/10 12 cm, and T 4 = 
T/10 4 K. Recalling that multicomponent effects are important 
for x p i > 0.1, then inserting this value into (|9} we can obtain 
the minimal metallicity Zm for which one-component models 
are applicable. By application of (JHJ in (|9} we find that 
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The latter equation can be rewritten in terms of basic stellar 
parameters using a canonical estimate of the wind terminal 
velocity v 2 ^ — a/(l — a)v 2 sc (see Castor et al. 119751 here- 
after CAK), where v csc is the escape velocity. Thus, the lowest 
metallicity for which one-component wind models can be ap- 
plied scales as 



Z, 
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In any case, the general result from Ea. il 1> is that multicompo- 
nent effects are likely to play a role only in very low metallicity 
massive stars. 

3. Model description 

To test this scaling analysis we calculated multicomponent 



models of stars initially modelled by Kudritzki (2002) and 
looked for multicomponent effects. In order to obtain cor- 
rect models for massive stars at extremely low metallicity and 



to improve the model convergence we used slightly different 
models from KKII. These changes are not important for the fi- 
nal results since these models retain basics characteristics of 
the KKII models. 

We assume that a stellar wind consists of three compo- 
nents, namely absorbing driving ions, passive nonabsorbing 
ions (hydrogen and helium, which mostly contribute to the bulk 
wind density) and electrons. For absorbing and nonabsorb- 
ing ions we solve the continuity equations, momentum equa- 
tions, and energy equations. The most important terms in the 
momentum equations are the radiative and frictional acceler- 
ation, while the most important terms in the energy equation 
are the heat exchange between wind components and the GO 
heating. Radiative force is calculated using the Sobolev ap- 
proximation (CAK) with finite disk correction factor (Friend 
& Abbott fT956l Pauldrach et al. 1986|. Contrary to KKII, for 
the electrons we solve the energy equation only. The radiative 
heating/cooling term in this equation is calculated using ther- 
mal balance of electrons (Kubat et al. [T999l Stellar fluxes at 
the lower boundary of the wind were determined using models 
calculated by a code of Kubat (2002). Electron density and ve- 
locity are calculated from the conditions of electrical neutrality 
and zero current. The charge of the nonabsorbing component 
is calculated using ionization balance of hydrogen and helium. 
We included helium to obtain the correct electron density. This 
proved to be crucially important in getting correct results for 
stars near the Eddington limit. Finally, we included new line- 
force multipliers introduced by Kudritzki ( 2002 ) both into the 
solved equations (with inclusion of the critical point condition) 
and into the linearization matrix (see Krticka 2002 ). 

The output of these models is the hydrodynamic structure 
of a three-component radiatively driven stellar wind, i.e. den- 
sity, velocity, temperature and charge of all components. For a 
more detailed description of these models see KKII. 

Originally, for the calculation of the radiative force CAK 
introduced a parameter related to the thermal velocity, calcu- 
lated at the stellar effective temperature. However, the radiative 
force itself does not depend on the thermal velocity. KKII as- 
sumed an artificial dependency of this thermal velocity param- 
eter on the local temperature instead of on the stellar effective 
temperature. To come closer to models calculated by Kudritzki 
( 2002 ) and to obtain a correct value of radiative force we as- 
sumed that the thermal velocity parameter used for the calcula- 
tion of radiative force is constant. 

The temperature dependence of the radiative force has, 
however, an important impact on the calculated models. If the 
radiative force does not depend significantly on the wind tem- 
perature (as was consistently assumed in this paper) then the 
velocity difference between wind components can reach the 
sound speed and runaway instability occurs. On the other hand, 
if the radiative force is decreasing with increasing temperature 
(as was assumed by KKII) then only heating occurs, and is 
accompanied by the lowering of the terminal velocity, or by 
the backfall of hydrogen for stars with extremely low-density 
winds. 
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Fig. 2. The calculated model of a star with logL/L© = 6.57, 
T cff = 50 000 K and Z/Z Q = 0.001. The stellar wind temper- 
ature structure is modified by frictional and GO heating. The 
temperature of absorbing ions is slightly higher than the tem- 
perature of other wind components. Note that velocities and 
densities of nonabsorbing component and electrons are nearly 
the same. 



4. Calculated models 

Generally, calculated models where multicomponent effects 
are important can be divided into three groups according to the 
influence of these effects. Each of this model groups will be 
discussed separately. 

4. 1. Modified temperature structure 

For stellar winds with only mild maximal velocity difference 
(f r i — v rp ) I \j^f- ~ . 1 the only effect of the multicomponent 
flow is a slightly modified temperature structure compared to 
one-component models. Because the radiative force for models 
in this paper does not depend on the temperature, the outflow 
velocity is the same as for the one-component models. This 
is a difference compared to KKII where for this type of stars 
lowering of the outflow velocity occurs due to the lowering of 
the radiative acceleration. 



An example of such models is given in Fig|2] In this case 
the temperature structure is affected by the frictional and GO 
heating. However, the increase of the wind temperature is 
lower, compared to models of main-sequence stars with sim- 
ilar velocity difference, mainly due to the higher wind density 
and consequently larger radiative cooling. The wind velocity is 
not changed compared to the one-component case due to the 
neglected temperature dependence of the radiative force. 

4.2. Runaway instability in the outer region 

For stars with lower metallicity the wind density of absorb- 
ing component is lower and, thus, the velocity difference be- 
tween wind components is higher in order to maintain com- 
mon flow (see Eq.(|6}). The lower the density, the higher the 
frictional heating. However, this picture changes for extremely 
low-density winds when the velocity difference is comparable 

to the sound speed, (ty; — v rp ) / ' \J^- ~ 1. The wind is not 
stable in this case for the ionic Abbott waves any more (see 
Owocki & Puis 120021 Krticka & Kubat 12002b and the runaway 
instability occurs. This case did not occur in KKII because the 
temperature dependence of radiative force prevented the veloc- 
ity difference to attain higher values comparable to the sound 
speed. 

One of the models where such instability occurs is given in 
FigH Note that the outer model boundary is selected in such 
a way that the possible region of runaway instability is down- 
stream. 

4.3. Possible fallback of nonabsorbing component 

For stars with even lower density, decoupling occurs in the 
region where the velocity of the nonabsorbing component is 
lower than the escape velocity. This causes either the fallback 
of hydrogen and helium onto the stellar surface or possible cre- 
ation of some kind of clouds above the stellar surface (simi- 
lar to that suggested by Porter & Skouza 1999 1. Note that, in 
this case, time-dependent modeling is necessary to calculate 
the proper structure of a multicomponent stellar wind. 

These conclusions are demonstrated in Fig|4] The absorb- 
ing component dynamically decouples from the nonabsorbing 
component, accompanied even by decoupling of the tempera- 
tures. 

4.4. Formula for the maximal velocity difference 

To test the reliability of our approximate formula for the max- 
imal velocity difference (jfji we compared the approximate ve- 
locity differences with model ones for stars with negligible 
multicomponent effects. The results are given in Fig[3] Clearly, 
formula (|6jl gives faithful predictions for the maximal drift ve- 
locity between absorbing and nonabsorbing wind components. 

5. Conclusions 

For most of the stellar wind models (moderate metallicity) cal- 
culated by Kudritzki (2002i we did not find any significant 
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Fig. 3. Calculated model of a star with logL/L© = 6.42, 
T cff = 40 000 K and Z/Z Q = 0.001. In the outer parts of the 
wind the velocity difference is comparable to the sound speed, 

(v r [ — v rp )l \J^~ ~ 1 an d runaway instability occurs. Thus, 
the outer model boundary is below this area. The GO heating 
has a negative sign near the stellar surface and positive in the 
outer wind regions. Due to this variation the temperature of ab- 
sorbing ions is lower than the temperature of other components. 

multicomponent effects, our calculations being in good agree- 
ment with his. Thus, radiatively driven winds of these stars can 
be adequately described by one-component models as done by 
Kudritzki. However, for some low metallicity stars we found 
that the multicomponent wind nature has a large effect and 
must be included to describe stellar winds of such stars. The 
importance of multicomponent effects increases with decreas- 
ing metallicity and for any given star we are able to find a 
metallicity value below which the multicomponent effects sig- 
nificantly influence the wind structure. We found several con- 
sequences of multicomponent effects, ranging from modified 
temperature structure to runaway instability and possible fall- 
back of the nonabsorbing component. Compared to results ob- 
tained by KKII, the frictional and GO heating are not so effec- 
tive mainly due to the larger radiative cooling of higher density 
winds. Note that the changed wind structure influences the stel- 
lar radiative flux and may be specially important in the UV or 
X-ray region. Finally, due to their different terminal velocities 
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Fig. 4. Calculated model of a star with logi/L© = 6.30, 
T cff = 40 000 K and Z/Z Q = 0.001. For this star the absorb- 
ing component is not able to accelerate the nonabsorbing com- 
ponent and the nonabsorbing component either form clouds 
around the star or falls back to the stellar surface. The increase 
of temperature is caused by the frictional and GO heating. Note 
that velocities and densities of nonabsorbing component and 
electrons are nearly the same. 
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Fig. 5. Comparison of the maximal velocity difference calcu- 
lated with approximate formula l|6) and model ones. 
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these stars cannot be used for stellar distance measurements 
using the wind momentum-luminosity relation (cf. Kudritzki 
et al.[T333}. 

The validity of models presented here is influenced by some 
simplifications used. Although they do not influence the ba- 
sic results, the real picture of multicomponent flows may be 
slightly different. Firstly, the radiative force is influenced by 
the temperature due to the temperature dependence of ioniza- 
tion and excitation. Although in the case of line-driven winds 
the occupation numbers are mainly given by the radiative pro- 
cesses, the temperature may be important especially in the case 
when large frictional heating occurs. For precise calculation 
of the radiative force it would be necessary to solve NLTE 
rate equations in the stellar wind. Secondly, the radiative heat- 
ing/cooling term may be influenced by complex ionization and 
recombination processes. Finally, for extremely-low density 
stellar winds there may be an insufficient amount of metallic 
optically thick lines and H, He lines may become important. 
We plan to adress these interesting issues using our NLTE wind 
code. First results obtained using this code for normal stars are 
promising and will be published elsewhere (Krticka & Kubat 
2003). 
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